1 




NASA CONTRACTOR REPORT 177340 


i 



i 


STUDY FOR PREDICTION OF ROTOR/ 
WAKE/FUSELAGE INTERFERENCE 
PART II: PROGRAM USERS GUIDE 


L 


(NASA-Cfi-177340-Vol«2) STOnI P 0a PREDICTION N85-?2347 

OF BOTOR/HAKB/FUSELAGB INTBBFBESNCE . Piai 
2: PBOGEAH USEES GUIDE Final Seport, 1 

Jon. 1980 - 1 Nov. 1983 (Analytical Hethods, Unclas 

Inc., Eediond , Nash.) 94 f HC A05/HF A0 1 G3/u1 20170 


! D. R. Clark 

| B. Maskew 

i 




/ 




NASA CONTRACTOR REPORT 177340 



STUDY FOR PREDICTION OF ROTOR/WAKE 

FUSELAGE INTERFERENCE 

PART II: PROGRAM USERS GUIDE 


D. R. Clark 
B. Maskew 

Analytical Methods, Inc. 
2047 - 152nd Avenue, N.E. 
Redmond, WA 98052 


Prepared for 

Ames Research Center 

under Contract NAS2-10620 


NASA 

National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett Field. California 94035 




TABLE OF CONTENTS 


Section 

LIST OF FIGURES 


Paqe-Mfl- 

iii 


1.0 INTRODUCTION 

2.0 GENERAL PROGRAM ARRANGEMENT * 

3.0 PANEL MODEL DEFINITION (Card Sets 9 through 16) 

3.1 Body Panelling 

3.1.1 Conventional Input 

3.1.2 Bodies of Revolution 

3.1.3 Lifting Surfaces — Type-3 Patches . . . . 

3.1.4 Rotors— Type-4 Patches 


4.0 WAKE INPUT 

4.1 Wake Grid Planes; Card Sets 17 and 18 .... 

4.2 Separation Line Specification; Card Sets 19 

through 23 

4.2.1 Wake Definition for Type-1 Wakes . . . . 

4.2.2 Wake Definition for Type-4 Wakes on 

Bodies ... 

4.2.3 Wake Definition for Type-4 Wakes on 

Rotors 

5.0 ROTOR BLADE ELEMENT MODEL INPUT DESCRIPTION . . . . 

5.1 Rotor Patch Identifiers (Cards R1 and R2) . . • 

5.2 Output Print Controls . - 

5.3 Iteration Controls (Card 14) 

5.3.1 Blade Flapping .... 

5.3.2 Rotor Lift 

5.3.3 Rotor Moments 

5.4 The Blade Element Model Geometry (Card Sets R6, 

R7, R8, R9 and R10) 


30 

30 

33 

42 

45 

51 

51 

51 

52 

52 

53 
53 

53 





TABLE OF CONTENTS (CONCLUDED) 


Section 


5.5 Rotor Performance and Control (Cards Rll, R12, 
R13 and R14) 


5.5.1 Rotor Speed and Flapping (Card Rll) . . 55 

5.5.2 Rotor Blade Cyclic Control (Card R12) . . 55 

5.5.3 Rotor Loads and Moments (Cards R13 and 

R14) 56 

5.5.4 Airfoil Data Sets (Card R15, Card Sets 

R16) 56 


6.0 INPUT DATA DECK BLOCKING AND VARIABLE LIST 57 


6.1 Input Summary 

6.2 Input Variable List 

6.3 Rotor Input Data Deck Description 


7.0 REFERENCES 73 


APPENDIX: Sample Case 


(a) Input 

(b) Output (Abridged) 






J> 


ft JF - £r JT ’‘ Jgfr 


-M ’M 


LIST OF FIGURES 


Title 


Page Ho . 


General Program Schematic 

Rotor Calculation Schematic 

Conventional Body Input 

(a) Patch-Section-Point Breakdown 

(b) Typical Input Data Set 

Typical Wing Panelling 

(a) Point Input and Tip Closure Schematic . . * 

(b) Typical Wing Input 

Sample— Body of Revolution— Nacelle 

(a) Patch Section Points Schematic ...... 

(b) Nacelle Input List 

Simple Elliptical Fairing 

(a) Basic Set-up 

(b) Input Data Listing 

Sample Lifting Surfaces 

(a) Simple Tail — Schematic 

(b) Simple Tail — Input Details . . 

(c) More Complicated Tail — Schematic 

(d) More Complicated Tail Input Details . . . . 

Main Rotor Modelling 

(a) Schematic 

(b) Input Details 

Tail Rotor Modelling 

(a) Schematic 

(b) Input Details 

Wake-Grid-Plane Definition 


3 

4 


6 

7 


9 

10 


14 

15 


16 

17 


20 

21 

22 

23 


25 

26 


28 

29 


(a) Simple Isolated Rotor 

(b) Embedded Body . . . . 


31 

32 


:y w ; f •' v T* mb# - v '* ^ :*£$ r ;'-:.'*w*«w jSv^ ' wk T*m*£ • 


Pig- No. 


LIST OF FIGURES (CONCLUDED) 


Title 

Type-1 Wake-Wing Separation Lines 


Page No. 


Wake Trajectory Definition 


Wake Shedding Schematic — Multi-line Input 


Multiple Patch Wakes 


(a) Type-4 Wake on Nacelle 

(b) Type-4 Separated Wake on Block Model of 

Rotor Head ..... 


Rotor Wake Specification — Simple Isolated 
Rotor 

Rotor Wake Input Schematic with Fuselage 
Present ...... .... 

Rotor Blade Element Model — Schematic . . 


1.0 INTRODUCTION 


Because the body/rotor program is a direct development of 
Program VSAERO (Ref* l)r no attempt will be made here to describe 
in detail the set-up and operation of these features of the 
present analysis which are in common with the original program. 
For the description of the non-rotor components, the development 
of the panelling schemes and the discussion of the underlying 
theory, the reader is referred to the original reports. The 
discussion here will place the rotor elements in the basic frame- 
work provided by VSAERO and will concentrate on a description of 
how rotary wing aircraft may be built up and placed in space 
using the flexibility inherent in the program geometry routines. 

The body/rotor code outlined here is directly related to the 
1000-panel version of Program VSAERO delivered to the NASA Ames 
Research Center and installed on the CRAY computer. The computa- 
tion times are consistent with those published in the original 
Program User's Guide for cases with complicated wakes if allow- 
ance is made for the rotor calculation by multiplying by a factor 
of roughly 1.3. 
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2.0 GENERAL PROGRAM ARRANGEMENT 

The rotor analysis is contained in a subroutine of program 
VSAERO and is called whenever a component is identified as 
belonging to a type 4-patch. Since the rotor blade element 
analysis responds to changes in inflow velocity/ the rotor calcu- 
lation is placed within the wake relaxation loop structure. This 
is highlighted in the solution block diagram in Figure 1. A 
typical solution would proceed as follows. 

A. From input data, program assembles body and rotor panel 
model. 

B. Program forms surface panel influence coefficients. 

C. Initial wake shape is constructed. 

D. On first pass, program sets up blade element model, and 
assuming uniform inflow, calculates rotor loads and mo- 
ments and (if requested) trims rotor. On subsequent 
iterations, inflow calculated on previous cycle is used. 
From calculated loadings, program evaluates rotor disc 
flow through velocities (local momentum balance) and 
doublet strength (to feed the rotor wake) and passes 
these back to the panel model. 

E. Solve for body panel singularity strengths. 

F. Calculate new wake shape and determine new inflow veloc- 
ities at rotor panel centers. 

G. Return to D to complete wake relaxation cycle. 

H. If requested, carry out viscous/potential flow iteration 
re-entering at D. 

I. Termination. 


Embedded within the rotor performance module are nested 
loops which control the blade flapping, rotor thrust (with col- 
lective pitch change) and rotor moments (with cyclic pitch 
change). These are outlined in the rotor calculation block 
diagram in Figure 2. At each azimuth location, the section loads 
are calculated for every station out along the blade. These are 
integrated radially to form the azimuthal totals. The blade is 
then moved to the next azimuthal location with flapping motions 
in response to any out of balance at the first position added. 
The blade is cycled around the azimuth until blade flapping has 
stabilised. The thrust is checked against the desired level and 
the collective pitch adjusted. Once the required thrust has been 
achieved, rotor moment trim is checked and adjusted if required. 
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GENERAL CONDITIONS 
SURFACE GEOMETRY 
HAKE GEOMETRY 



AND MOMENTS I THROUGH THREE VALUES 

i • TWO DIRECTIONS 

OFF-BODY VELOCITY SURVEY . UNEQUAL INTERVALS) 














3.0 PANEL MODEL DEFINITION (Card Sets 9 through 16) 


3.1 Body Panelling 

3.1.1 Conventi onal Input 

Body panelling, including any lifting components wh 
present, follows exactly the format set down in the Progr 
VSAERO User's Guide, Ref. 1. As a result, the body can be 
up of any convenient combination of panels, patches and com- 
ponents. In the VSAERO hierarchy, the panel is the smallest 
element controllable by the user. A patch is a collection of 
panels in a regular array and a component is an identified group 
of patches taken together for the evaluation of force and moment 
totals. Components may be further grouped into assemblies. 

The body shape may be entered in several ways, depending on 
its complexity and on the way in which the shape data is avail- 
able. Conventionally, the body is defined with a set of loft 
lines which are generally cut at constant body locations, station 
(x), buttline (y) and waterline (z). Figure 3 illustrates how a 
typical body may be first broken into patches, Figure 3 (a), and 
then input. Figure 3(b). The patch allocation is made so that 
the different regions of the body may be represented by panelling 
of the appropriate density, high in regions of particular in- 
terest, low in other areas bearing in mind that within a patch, 
the panels form a regular array with the same number of panels in 
each row and column. Consequently, patch boundaries almost al- 
ways occur where large changes in body cross section are present 
and in regions where, away from areas of interest, panel densi- 
ties are being reduced for reasons of economy. 

In the case of the body used in this illustration, the shape 
was defined with a series of station (constant x) cuts, fol- 
lowing the VSAERO User Document, Ref. 1, each section is defined 
by a series of points located in a local axis system in the 
global coordinate system. The data set defining each section 
contains a header card which contains the origin of the local 
axis system in the global coordinate system, the scale and orien- 
tation of the section (section may be scaled and/or rotated to 
ease input), and indicator cards which alert the computer to the 
way in which the data is being input and which indicate whether 
the section is internal to or closes a patch and provide informa- 
tion on how the patch is to be divided up. The header card is 
then followed by the string of points defining the section, 
together with node cards which alert the computer to the way in 
which the surface is subdivided, to any changes in surface curva- 
ture at junctions, and to the termination of the string. 

Program VSAERO offers the user great flexibility in the ways 
in which the shape may be input. Individual defining sections 
may be input in either the x (with y and z), y (with x and z), or 
z (with x and y) planes or in generalized x, y and z coordinates. 
Any other previously defined section may be automatically copied 



Patch-Section-Point Breakdown. 
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and a whole family of automatically defined airfoil, elliptical 
or polar coordinate sections are available. 



The data sets for successive sections are stacked to form 
the input data set for a patch and the way in which they are 
stacked is of particular importance. Sections should be stacked 
in order and in a direction consistent with the way in which the 
points are input along each section. In the example shown in 
Figure 3, the points on each section were input up along the 
profile on the start ard (positive y) side toward the top center- 
line and although the input points do not necessarily dictate 
panelling directly, the direction of input determines the order 
of the panel, with the first section input making up side 1 of a 
patch, and the panels numbered in the direction of the input 
points. Since program VSAERO requires that panels and patches 
have an anticlockwise corner point order when viewed from outside 
with the bottom to top input scheme used in Figure 3, the sec- 
tions must be stacked from front to rear. If the points had been 
entered from top to bottom, tnen the section input order v/ould 
have been reversed and the patch would have been defined with 
sections input first at the rearmost edge of the patch and 
working forward. 

Input sections need not have the same number of defining 
points. However, when instructing the program on how the section 
should be subdivided, it should be remembered that a patch is a 
regular array and all columns must have the same number of 
panels. This is determined by the user properly setting the 
appropriate dividing instructions (the node cards) in the section 
input string. 

Wing sections are conventionally input as shown in Figure 4. 
Using the same format as outlined for bodies above and starting 
with the lower tr ailing-edge point, the surface is input point by 
point working forward toward the leading edge, around the leading 
edge and aft over the upper surface. As with the body input 
described above, node cards are inserted to delineate regions of 
differing panel density or changes in surface curvature. For 
instance, in the example shown, the panelling is required to be 
more dense close to the leading edge. This is achieved by in- 
serting a node card with a value of N0DEC*1 at the leading edge 
(indicating the end of a region but maintaining surface curva- 
ture) and using values of the distribution controls, INTC-2, at 
the end of the lower surface and INTC-1 at the end of the upper 
surface. Following this procedure, sections should be stacked 
from root to tip if only the starboard (positive y) side of the 
vehicle is being modelled. 

Although a panel is the basic element in modeling the sur- 
face and the solution proceeds with the strength of the singu- 
larity used to represent the panel as the unknown, the patch is 
the most conveniently manipulated unit. Patches may be construc- 
ted and moved into place to represent the configuration in any 
way convenient to the user. They may be ordered in any way 
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Point Input and Tip Closure Schematic. 
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without significantly affecting the solution and allow changes to 
a configuration to be made by removing a feature and simply 
plugging in one or more replacement patches to remodel the 
vehicle. 

Summarising, the input list for a typical surface patch 
would be as outlined below. The names used for the data items 
are those employed in the Program VSAERO User's Manual. The 
terms "chordwise” and "spanwise" can be considered to represent 
directions along the input section and along the patch length 
respectively. For completeness, they are also defined here. 


Patch Card :- CARD IQ . 

IDENT Patch type 1-Wing 

2- Body 

3- Lifting Surface (Neumann) 

4- Rotor 

KOMP Defines component and assembly to which 

KASS patch belongs 

PNAME Patch title 


Section Cards :- (Stacked in spanwise direction) CARD 11 
STX, STY, STZ Location of origin of section input points 

SCALE, ALPHA, BETA Scaling factor and rotation angles 
INMODE Type of section input 



NODES 

Nodes to signal changes in spanwise panel 
distribution, surface curvature, end of 


patch, etc. 


NPS 

INTS 


Number of spanwise panels in the interval 


Way in which spanwise 


section is divided 


Defining I nput Points :- CARD SET 12 . 

BX, BY, BZ Coordinates of section points. Form depends 

on INMODE used. 


\ 


* 


11 


n 


CARD 14 


NODEC 


Nodes to signal changes in chordwise distri- 
bution, surface curvature, end of section, 
etc. 


INTC 


Number of chordwise panels to be generated in 
in the interval 

Way in which chordwise section is to be 
divided 


The values of the nodes used to delineate changes in either 
chordwise (NODEC) or spanwise (NODES) directions are the same. 
In both cases, values of 1 and 2 represent the end of a region 
within a patch where in the first case, surface curvature is 
continuous into the next region and in the second case is discon- 
tinuous. A value of 3 indicates the end of a section, chordwise, 
or a patch, spanwise, while values of 4 and 5, used only on the 
last sections of patches indicate, respectively, the final sec- 
tions of components and of the whole configuration. 

A similar arrangement with matching values is used to indi- 
cate the way in which the surface is to be divided in the 
chordwise and spanwise directions. The parameters, INTC and 
INTS, may be values of 0 through 3 depending on whether the 
points are to be: 

Value 

Closely spaced at the beginning and end of the region 0 
Closely spaced at the beginning of the region 1 
Closely spaced at the end of the region 2 
Equally spaced 3 

When an unequal discribution is called for, the distance 
along t he surface across the interval is broken up using a cosine 
form. A detailed description of the procedure is given in the 
Program VSAERO Manual, Ref. 1. 

Of course, the panelling could be built directly upon the 
input points and sections with no subdivision or interpolation. 
In this case, the input points on successive sections are simply 
connected together to form the panelling. This option is selec- 
table by setting NPS or NPC equal to zero* Any values other than 
zero determine the number of panels within the particular span- 
wise or chordwise interval. 
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In the discussion above* the standard method of inputting 
configurations has been outlined. Program VSAERO offers many 
other options. Two of them* which are particularly useful in 
entering helicopter shapes, are outlined below. 


3.1.2 Bodies of Revolution 

There are many applications in the modelling of typical 
helicopter shapes that do not require the detail provided by ( 

section by section input. Many configuration elements* in fact* 
may be simply defined as bodies of revolution. Examples are 1 

engine nacelles* tail boom sections* rotor head fairings and 
external stores. Figures 5 and 6 illustrate how two of these, a 
simple engine nacelle and a rotor head fairing (or radome) may be 
generated. V 

For the first example* Figure 5* the engine nacelle* the 
unit is made up of one section defined in the y * 0 plane by 
inputting values of x and z to outline the profile and then 
rotated through 3600 to form the nacelle. It can be modelled as 
one patch with four separate regions. They ares the inlet face* 
the inner surface of the inlet* the outside surface of the na- 
celle* and* finally* the exhaust plane. 

As with normal input, the card set describing the nacelle is 
preceded by a patch card. Generally* only one section card is 
required. This locates the input section in the global coordi- 
nate system and specifies its orientation and scale relative to 
that frame. The parameter, INMODE* describes how the section 
points are to be input; INM0DE»2* for both the examples in 
Figures 5 and 6 signifies points being entered in the y ■ 0 
plane. If a flow-through nacelle were being modelled* INM0DE»5 
could have been used to specify a standard NACA 4-digit cross 
section. 

The body of rotation option is selected by setting the 
spanwise node card* NODES* to a negative value. For the ex- 
amples* simple polar symmetric bodies are used so NODES is -3 in 
both cases. However* more complex bodies may be built up by 
combining sections to vary panel density* and in that case* 
values of -1 or -2 would be appropriate at the intermediate 
spanwise sections depending on whether the surface curvature was 
continuous or discontinuous across the node. The final section* 
completing the rotation* would* of course, have NODES— 3 signi- 
fying the closing of the patch (or -4 or -5 if this were the last 
patch on the component on the body) . 

Definition of the cross-section points (defining the local 
chordal shape) follows the same rules for defining conventional 
cross sections outlined above. 
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Figure 6. Simple Elliptical 
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(b) Input Data Listing. 
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The input set is completed by specifying the final and 
initial angles, a, to define the rotation arc on CARD 15. For a 
body such as the nacelle, remote from the system plane of sym- 
metry, a full 3600 of rotation is used. However, in the same 
applications where the flow is symmetrical, where no yawing is 
involved and where the body is on the plane of symmetry, then 
only the positive side of the body of revolution is defined and 
is 1800 . 

The second example used, the ellipsoidal fairing of Figure 
6, illustrates another feature of Program VSAERO; that is, the 
ability to construct elements of the configuration as separate 
components in the most convenient orientation and then rotate 
them or move them to the required location. The positioning and 
rotation are controlled by entries on a special component card 
which must precede a patch or a group of patches identified as 
belonging to a separate component. This card precedes the first 
patch identifying card. 

In the simple example given the component is made up of only 
one patch. As with the nacelle, it is defined in the y * 0 plane 
by inputting the values of x and z outlining the cross section. 
This outline is then rotated about the section 'X' axis to form 
the body of revolution. 

Once the shape is formed in the section coordinate system, 
three further processes are available at the component level 
before the final shape is set in the global frame. At this 
level, all the patches within the component may be scaled, either 
increased or reduced in size, rotated or translated. The changes 
are applied in that order. Two rotations are available; a simple 
one, where the body is rotated about the component 'Y' axis and 
rotation about a general axis that is defined by two input 
points. The simple rotation is the default mode. The user- 
defined axis rotation is entered by specifying a negative value 
for the scale parameter and then inserting an extra data card 
which contains the rotation axis defining points. For the simple 
fairing case shown in Figure 6 the default is used and the body 
is pitched up 90o to lie with its long axis in the horizontal 
plane. 


3.1.“* Lifting Surf aces— Tvpe-3 Patches 

A further option offered by Program VSAERO which is 
useful in modelling components remote from a region of interest 
is the ability to use a lifting-surface representation rather 
than a full surface singularity model, and enforcing the Neumann 
rather than the Dirichlet boundary condition. A typical applica- 
tion of this feature would be in modelling wings or other lifting 
or control surfaces that do not have a direct effect on the com- 
ponents that were the focus of the study or in regions where no 
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aerodynamic rigor is sacrificed by not including the effects of 
the thickness of the components. Figure 7 shows examples of this 
where simple type-3 patches are used to model typical 
vertical/horizontal tail assemblies. 

Depending on the level of detail required, this assembly 
could be created by as few as two defining points on each of 
three sections if simple, flat surfaces are desired. The example 
in Figure 7(a) and (b) shows this input. The group of input 
cards required is preceded by a patch card (and a component card 
if a separate component is called for) calling out a type-3 
patch, input continues with the section card defining the lower 
root and then two cards defining the trailing- and leading-edge 
points followed by a chordwise node card. Note that on the 
chordwise node card the program is instructed to break the chord 
down into five panels, NPC-5, distributed so that they are denser 
at the leading edge, INT02. 

The second section card, at the joint between vertical and 
horizontal planes, with a spanwise node value, NODES, of two 
signifying a discontinuity in surface curvature, contains the 
information to divide the vertical section in the spanwise direc- 
tion into four rows of panels, NPS»4, spaced equally, 
INTS»3. Input is complete with the section defining the horizon- 
tal plane. In this case there are five panels spanwise, NPS»5, 
divided so that they are spaced densely at root and tip, INTS*0. 
Since this is the last section to be input on this patch, the 
spanwise node, NODES is set equal to three. The panelling in the 

chordwise direction is the same on the horizontal plane as in the 
vertical. 

In the second example. Figure 7(c) and (d), cambered sec- 
tions are used on both the vertical and horizontal sections and, 
consequently, a small transition piece is required to go from the 
vertical to horizontal planes. This is formed quite simply by 
allowing the program to correct the sections defining the top of 
the vertical panel and the root of the horizontal panel and 
generate a row of panels to fill the gap. To do this, four 
sections are required, The first section, as before, defines the 
root of the vertical panel this time using INMGDE«1, points input 
as X and Y values, to define the section camber line. The sec- 
tion was defined relative to the quarter chord point; tha values 
of STX, STY, and STZ locate the section in the global frame. In 
the example, the camber line distribution is assumed constant up 
the panel and as a consequence, the top section is input by 
simply relocating the root section with the appropriate values of 
INMODE a O, with the SCALE factor set to give the correct taper. 

The root section of the horizontal plane is again defined 
with input points, this time values of X and Z in a constant Y 
plane, INH0DE>2, taking care to place it in the correct position 
relative to the earlier input sections using the appropriate STX, 
STY, and STZ local origin values. Again the tip section (the 
camber is assumed constant) is input simply by copying and 
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scaling the toot section. Hete, however, the section is pitched 
down using the pitch angle control, ALF, to model an appropriate 
value of section twist. Since both surfaces and the connecting 
fairing are being modelled as one patch, the chordwise distribu- 
tion of subdivided panels must be conserved. This does not 
require, as in this example, that the same number of points be 
input provided that the chordwise node cards for each section 
have a non-zero value for NPC, the number of panels chordwise, 
and that this is repeated at each spanwise station. 


I 

I 


3.1.4 RQtQcs-Type -4 Patches 

In the body/rotor analysis, two models of the rotor are 
constructed. They are the detailed blade element model which is 
to be discussed later and the panel model which provides the 
coupling between the rotor and fuselage flow fields. The panel 
model of the rotor is constructed in much the same way as the 
bodies of rotation discussed above with the added .act that the 
chordwise and spanwise (in the rotor framework radial and azi- 
muthal) panel breakdown forms the basis for the radial spacing of 
blade stations and the azimuthal increments in the blade element 
calculation. 

The rotor is modelled using a disc-like array of panels 
generated in exactly the same way as the bodies of revolution. 
The panels are represented, as are all panels in program VSAERO, 
by a combination of source and doublet singularities. In a type- 
4 patch both source and doublet elements of the singularity are 
specified (known as a result of the blade element calculation) 
and are passed over, within the program, from the blade element 
to the panel calculation. Part I of this report contains an ex- 
planation of the form of the singularities. Each rotor disc must 
form a separate patch preceded by the appropriate pacch card and 
it is recommended that they also be identified as a separate 
component for ease of manipulation and the separate accumulation 
of loads. Figure 8(a) and (b) shows how a typical main rotor may 
be generated. 

The input for this example. Figure 8(b), begins with a patch 
card with IDENT set to four (4) to signify a rotor. This is 
followed by the section card, placing the rotor center in the 
global frame with the appropriate values of STX, STY and STZ and 
announcing that a body of revolution is to be generated with 
HODES«- 3, -4 or -5. Using this approach, INMODE is most commonly 
two (2) and the line defining the rotor disc is input in the y - 
0.0 plane with values of x and z and rotated about the x-axis 
(default). It is at this stage that blade coning can be intro- 
duced, as in the example. With this simple form only two points 
are required, but, of course, any higher-order description may be 
used to generate curved rotor disc surfaces. In the example, the 
point closest to the axis forms the innermost defining station 
for the aerodynamic parts of the blade and the outer point the 
tip radius. Normally the blade radius is broken down 
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Figure 8. Main Rotor Modelling. 









Input Details. 


automatically specifying/ as in the example/ the inner and outer 
points and the number and distribution of panels to be formed on 
the chordwise (radial) node card. INTC-2 is recommended for the 
radial breakdown since this concentrates the panels at the outer 
end of the radius. The number of panels is dictated by the 
detail required. Cases have been successfully executed with a3 
few and 3 panels radially and the maximum set by program capacity 
is 20. Experience with conventional rotor performance programs 
has shown that fifteen segments radially is more than adequate 
for detailed blade studies. 

Since the disc is being formed by the body of revolution 
feature/ the number of azimuth stations is set by the NPS pa- 
rameter on the section card and the distribution should be uni- 
form with INTS set equal to three (3). The input is completed by 
including the rotation angle range of 3600 . 

With the disc panelled in the input location about the x- 
axis/ all that remains is to pitch it into the appropriate atti- 
tude. This is accomplished using the rotations available on the 
component card. The default rotation is about the component y- 
axis with a positive rotation being nose-up. In the example/ it 
was determined that the tip path plane should be 50 nose-down so 
the rotation required on the component card to place the disc in 
the correction orientation was +85o. The values of the component 
origin, CTX, etc. then locate the rotor relative to the rest of 
the configuration 

If a more involved rotation was required, say to apply some 
lateral tilt to a main rotor or to position a tail rotor, this 
can be done by activating the option to specify the axis of 
rotation. This is done by setting the scale parameter on the 
component card to a negative value and inserting on the following 
card two points specifying the rotation axis. The details for 
applying this technique are discussed in full in the VSAERO 
program guide and are summarised in Figure 9(a) and (b) for the 
case of a tail rotor. Here, as a further example of the program 
flexibility, the rotor is input with unit radius and the scale 
parameter on the component card is used to bring the radius up to 
the full scale value. A further point to note in the tail rotor 
model is that the rotation angle range has been changed to: 8 2 • 
450.0, * 90.0. This is required to line up the first panel in 
the* panel model of the rotor with the zero azimuth location in 
the blade element model. Since this is conventionally oyer the 
tail for the main rotor and along the aft pointing horizontal 
radial for a tail rotor, the tail rotor panel model before rota- 
tion must start in the horizontal position. The same effect 
could, of course, be achieved by using INMODE-3 or 4, carrying 
out the original line definition in the z - 0 plane and by 
returning to the « 2 - 360.0, 81 - 0-0 degrees rotation range. 
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4.0 WAKE INPUT 

4.1 Bake Grid Planes :- CARD SETS 17 and 18 

In program VSAERO the development of the wakes after they 
are shed is calculated in a series of planes perpendicular to the 
onset flow starting upstream of the first shedding location. The 
way in which these planes are set up is described in full in the 
VSAERO program guide and the recommendations made there for fixed 
wing aircraft and general shapes apply equally well to rotor/body 
problems. 

Figure 10(a) illustrates how the defining grid planes would 
be set up for a typical helicopter study. Note that the first 
grid plane corresponds to the leading edge of the disc. As a 
general guide for isolated rotors f the grid planes should be 
distributed over the disc to correspond to the azimuthal break- 
down of panels using a full cosine distribution (dense at 
leading- and trailing-edges, open in the center) with steps equal 
to half of the total rotor azimuthal increments. Downstream of 
the trailing edge of the disc, the first two disc diameters of 
distance could be broken into ten segments with half cosine 
spacing and then a further six diameters with four large segments 
to fill in the far-field effects. 

In regions where substantial flow distortion is expected, 
say around the front of a fuselage or if details of the passage 
of a wake over a wing leading edge are required, then regions of 
more closely set grid planes are required. Figure 10(b) is an 
example of how this may be achieved. 

As with the input of body sections, the wake stations are 
specified by inputting a series of values separated with node 
cards which indicate how intervals are to be divided up. The 
number of divisions may be specified using NPC; the distribution 
is specified as with body input using the parameter, INTO, and as 
before, individual stations may be input and used without further 
subdivision by setting NPC equal to zero. The last wake grid 
plane should be followed by a node card with NODE set equal to 
three (3) to terminate wake grid plane input. 


4.2 Separati on Line Specification !- CARD SETS 19 through 21 

In the version of program VSAERO described in this report, 
two types of wake are available. These are the type-1 wakes, 
springing from lifting bodies, and the type-4 wakes, enclosing 
regions of flow with energy states higher or lower than the 
surroundings. The way in which these wakes are described are 
generally similar, the attachment process is identical and they 
only differ in that for the type-4 wakes the velocity jump across 
the jet sheet must be identified in order that the program may 
assign the correct vortex strength to the sheet elements. 
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Simple Isolated Rotor. 


STATION 



(b) Embedded Body. 
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While the uses of the type-1 wakes are clear cut, on wing 
and other lifting-surface trailing edges, the use of type 4 wanes 

requires some additional explanation for the hell . copt akes 
tion several elements of helicopter configuration® shed wake..* 
wher'e chSng« in energy level ere involved Hcst obvioue is the 
from the powerplant. Less clear, but nonetheless impor 
tant is the wake from the rotor head assembly. Since, in this 
the^hape is normal ly not modelled in detail and the drag 

coefficient of the assembly is generally k: "5""* 

model the overall effect of the component hy =onstr noting a bul 
model and attaching a type-4 wake along its aft facing e g 
The rotor model also uses a type-4 wake. 


Each individually identified wake is treated in much the 

same way as a patch J handled in the body input and, as in the 

case of Y the patch, it must be preceded by a ies the 

CARD 19. This card contains the information which identifies the 

wake type, IDENTW-1 for a regular wake or -4 for a jet 
n Ad® 1/ separated base/rotor wake; indicates whether the wake is 
held^^^d*^!: isrelaxed, IFLEXW-1 or 0, respectively; and a 

descriptive title. 


4.2.1 nef in i tion for Tvpe-1 Wak fifi 


The method of attaching the wake along the separation line 

is explained in detail in the VSAERO Program User * G ?if e ; r iai- 
following notes should be considered a supplement to the ongi 
nal, more detailed presentation. 


Th* separation line is applied to th« ,“ f l T n the example 
that the local attached flow comes from the l*fi. “ the esample 
shown in Pigure 11 the separation is parallel to side 2 of the 
patch. The string of panels to which the wake sheet is attached 
is identified on Card 20 as follows. 


(KWPACH) (KWSIDE) (KWLINE) (KWPANl) (KWPAN2) (INPUT) (NODEWS) 


1 

(or 0) 


6 

(or 0) 


Default for the set 
of panels alongside 
KWSIDE 


Using the default, 0, for 
KWPANl, KWPAN2 in cases 
where the string of panels 
is the complete set parallel 
to side KWSIDE is recommended 


KWPACH is the patch to which the wake is attached and KWSIDE the 
p«?h side parallel to and in the same direction as the shedding 
line If the sheet is attached along the patch edge, KWLINE 
takes the default value of zero. However, if the s h® et . v t® r j 
attached along an internal panel edge, this would be identified 
by counting, in this case along side 1, until the shedding panel 


A*Vl 
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WAKE-SHEDDING PANELS 










A 


was reached. This number* then* defines KWLINE. KWPAN1 and 2 
define the spanwise extent of the shedding; again the default 
values of zero give full span, and numbering proceeds across the 
patch in the direction of the shedding line where only part-span 
shedding is present. Card 20 is equivalent to the section card. 
Card 11* of the body input. 

Because INPUT is 2 on card 20 above, this card must be 
followed by CARD SET 21/22 defining the geometry of streamwise 
wake-line, LINE 1. As with the input of body section data, the 
wakes may be prescribed in a number of ways. For the current 
example, INPUT-2 indicates that the program expects data to be 
entered in X and Z coordinates with a local origin at the shed- 
ding point. As with the body input, wake filaments may be des- 
cribed with any level of detail desired and the wake segmented 
with node cards separating the different elements. For the 
simple case used here, if an initially rectilinear wake was 
required, it should be enough to prescribe a point in the far 
field downstream and a node card. If the wake had to pass an 
obstruction, say another body or a flap, a more detailed path may 
be prescribed. Both examples are outlined in Figure 12 below. 
In both cases the X and Z coordinates are relative to an origin 
at the trailing edge. CARD SET 21/22 is equivalent to CARD SET 
12/14 of the body input. 

When the input is complete the wake is terminated with 
another CARD 20: 


(INPUT) (NODEWS) 
0 3 


The string of wake-shedding panels End of this wake 

has already been defined hut another wake 

must follow. The 
final wake must 
have a 5 here 

This copies LINE 1 to LINE 2 to 
complete the wake geometry des- 
cription. (Note that INPUT>0 
can be used here if a different 
wake line geometry is required, 
ir. which case CARD SET 21/22 must 
follow for LINE 2) 


/ 


i 


Note : In the case of type-1 wakes (IDENTW-1 on CARD 19) it is 
possible to turn the wake over and attach it along side 4 of the 
patch — in this case it would start at the tip and move inboard. 
CARD 20 would then be: 
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(WAKE 2 NOT SHOWN FOR CLARITY) 


I Figure 12. Wake Trajectory Definition. 


I 



(KWPACH) 

1 


(KWSIDE) (KWLINE) (KWPANl) 

4 0 1 

(or 0) 


ORIGINAL PAGE IS 
OP POOR QUALITY 

(KWPAN2) (INPUT) (NODEWS) 

7 2 0 

(or 0) 


No change 



In the initial example, the wake input at the first section was 
simply copied across the span, but in many cases this would be 
inappropriate. A more involved example is given below. 

Changing the wake line geomtry in the middle of a patch is 
achieved by simply breaking the string of wake shedding panels 
into sets. This is possible since there is an opportunity to 
define a streamwise wake line at the beginning of each string of 
wake-shedding panels and at the end of the last string. This is 
explained by Figure 13. The input for these wakes would be as 
follows 

IDENTW 

CARD 19 10 0 Wing Wake 

CARD 20 1201520 


CARD SET 21/22 SWPX(l), SWPZ(l) 

for Line 1 INPUT*2 Format 

SPWX (2) , SWPZ ( 2) 


(NODEWS) (NWP) (INTO 
3 10 3 


I 

2 0 6 6 0 0 

KWPANl, KWPAN2 INPUT -0 

Copies Line 1 

into Line 2 < 

2 0 7 12 2 

INPUT for 
Line 3 


CARD 20 


CARD 20 
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c 

SIDE VIEW 



— — — — LINE © 

L- LINE © 

4 . 

LINES © AND © 

3 " LINES © AND © 

(SECOND WAKE) 


Figure 13. Wake Shedding Schematic--Multi-line Input. 



CARD SET 21/22 
For Line 3 


SWPX(l), SWPZ(l) 
SWPX (2) , SWPZ (2) 


CARD 20 


(NODEWC) (NWP) (INTO 
3 10 3 


2 3 

NODEWS* 3 
terminates 
this wake 


CARD SET 21/22 
for Line 4 


SWPX(l), SWPZ(l) 
SWPX (2), SWPZ (2) 


(NODEWC) (NWP) (INTO 
3 10 3 


CARD 19 


Flap Wake 


CARD 20 


CARD SET 21/22 
for Line 5 


SWPX(l), SWPZ(l) 
SWPX (2) , SWPZ (2) 


CARD 20 


(NODEWC) (NWP) (INTO 
3 10 3 


INPUT-0 Copies Line NODEWS-5 

5 (i.e., the previous for final 

line) into Line 6 wake 


If the wake separation line passes over a number of patches then 
a separate string of wake-shedding panels must be specified for 
each patch. (Multiple strings of wake-shedding panels may still 
be specified within a patch as shown above). This is illustrated 
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with the example shown in Figure 14. Following the sketch, the 
multipatch input for this wake becomes: 


KWPACH KWSIDE KWLINE KWPAN1 KWPAN2 INPUT NODEWS 
CARD 20 12 0 0 0 20 

CARD SET 21/22 FOR WAKE LINE 1 

CARD 20 22 0 0 0 00 

(Copies wake line 1 into 
wake line 2) 

CARD 20 31 0 5 6 20 

CARD SET 21/22 FOR WAKE LINE 3 

CARD 20 32 0 0 0 20 


CARD SET 21/22 FOR WAKE LINE 4 

CARD 20 00 0 0 0 25 

End of Wake Input 


CARD SET 21/22 FOR WAKE LINE 5 

WARNING; Before leaving the discussion of type-1 wakes, it is 
important to note that the panelling on the upper and lower sides 
of the wake must match. This is required so that the correct 
doublet jump across the wake may be properly evaluated and shed 
into the wake columns. 
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PATCH 2 


ORIGINAL PAQB ft 
OF. POOR QUALITY 



Figure 14. Multiple Patch Wakes. 
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4.2.2 Wake Definition for Tvpe- 4 Wakes on Bodies 

Type-4 wakes are specified in the same manner as type-1 
wakes, and are fed by flow from the left when looking along the 
direction of the separation line. Figure 15 shows two examples 
of this type of wake. 

In the jet efflux case. Figure 15(a), the nacelle and base 
are all part of the same patch and the wake is attached along the 
aft-facing edge. This is parallel to side 2 of the patch on the 
17th row of panels along side 1 (see Figure 5 for details of 
nacelle panelling). Consequently, KWPATCH takes the nacelle 
patch number, KWSIDE»2 and KWLINE»17. Since the wake is attached 
across the full width of the patch, KWPAN1 and KWPAN2 are set to 
zero. 

In the case of the rotor head block model. Figure 15(b), the 
block is in three parts with separate patches for the front and 
rear faces, patches 1 and 3, and the main surface as patch 2. 

Here, the wake is attached along the edge, side 3, of patch 2. 

Consequently, KWSIDE*3, KWLINE-O (the default for the edge) and 
KWPAN1 and KWPAN2 are again zero, since the wake goes all around 
the edge of the patch. 

The wake geometry for type-4 wakes is defined in a manner 
identical to that required for type-1 patches on CARD SETS 21/22. 

Since type-4 patches involve regions of higher or lower 
energy and the wake strength is set to produce the appropriate 
internal flow, the inner and outer velocities (normalized with 
respect to a unit onset flow) must be specified. This is done on 
CARD SET 23 immediately following the wake definition. 

Input for normal type-4 wakes (i.e., those attached to 
normal wing or body patches) is completed by ensuring that the 
flow into the wake cavity, from the area of panels enclosed by 
the wake, matches the flow down the inside of the wake column. 

This is done by using the option available within the program to 
suspend the usual assumption of zero flow through the boundary i 

and by replacing it with appropriate normal velocity. This would 
be positive for an outflow. The special options available with 
CARD SET 8 provide this capability. 

The first entry on Card 3, NORSET, sets the number of < 

regions in which the transpiration velocity is to be changed. If 
this is non-zero, the program then expects to read a Card 8A for 
each region. Explained in detail in the VSAERO manual, Card 8A, 
identifies the patch and the rows and columns involved, and 
specifies the /elocity value. Examples are given of typical CARD 
SET 8/8A input for the two cases in Figure 15. 


t- 
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KWSIDE=2; KWLINE=17 
KWPAN1 AND KWPAN2 = 0 


CARD 8 

(NORSET 

NVORT 

NPSOM 

JETPAN) 

1 

0 

0 

0 


CARD 8A 

(NORPCH NORF 

1 18 


NORL NOCF NOCL VNORM) 
1 0 0 2.5 


DEFAULT FULL WIDTH (TYPICAL) 


Figure 15(a). Type-4 Wake on Nacelle. 




ORIGINAL PAGE » 
OF. POOR QUALITY 




CARD 8 : 

(DORSET 

1 

CARD 8A : 

(NORPCH 

3 


Figure 15(b) . 
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NVORT ETC.) 
0 


NORF 


0 




NORL 

0 


NOCF 

0 


NOCL 

0 


VNORM) 

0.0 


DEFAULT FULL PATCH 


TYPICAL 


Type-4 Separated Wake on Block Model of Rotor Head 


-nr 





CARD SET 8/ 8A is also used to identify sets of panels which 
are known to fall inside regions of higher/lower energy level. 
This is required if correction of the calculated values of the 
pressure coefficients to account for the altered dynamic head of 
the region is desired. As with the outflow option above, the 
number of sets of panels to be modified are identified with 
JETPAN on CARD 8 and the patch row and column information sup- 
plied on CARD 8A for each set. The incremental dynamic head is 
calculated with the wake sheet values of VIN and VOUT input on 
CARD 23. 


4.2.3 wake Definition for. Tvpe-4 Wakes o: i. Rotors 

The wakes used in the rotor calculation are a special 
form of the type-4 wakes discussed above and are invoked by the 
program automatically when a type-4 wake is attached to a type-4 
rotor patch. Wake attachment is similar to that described for 
other type-4 patches above, but extra care must be taken to 
completely enclose the wake volume. This is detailed in Figure 


i 

i 


L 


Considering how the rotor patch is generated by rotating the 
input section, the chordwise (in this case radial) direction is 
side 1 of the patch. The outer edge of the disc becomes the 
spanwise direction and is side 2 of the patch. Side 3 and side 4 
follow naturally to complete the definition. The card set for 
the wake is as follows. 


£ ■ 


{ 


(IDENTW) (IFLEXW) ( IDEFW) 

(WNAME) 

CARD 19 4 0 0 

ROTOR WAKE TIP 

Type-4 Wake 



Separation line 

Distorted wake 

defined by panel 

will be calculated 

string. Card 20 


must follow 


(KW PATCH) (KWSIDE) (KWLINE) (KWPAN1) (KWPAN2) (INPUT) (NODEW-1) 
CARD 20 1200020 

Tip circle Default 

attaches Default attaches 
wake along wake along the 

patch edge full length of edge 


4 




f 
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CARD SET 21/22: DEFINING WAKE FILAMENT AS IN TYPE-1 WAKES ABOVE. 


(KWPATCH) (KWSIDE) (KWLINE) (KWPAN1) (KWPAN2) (INPUT) (NODEWS) 
0 0 0 0 0 0 3 

terminates 

copies filament input tip wake 

on previous CARD 20 to 
close wake column 


(IDENTW IFLEXW IDEFW) 

CARD 19 4 0 0 


(WNAME) 

Rotor Wake Root 


(KWPATCH) (KWSIDE) (KWLINE' (KWPANl) (KWPAN2) (INPUT) (NODEWS) 
CARD 20 1400020 


CARD SET 21/22: TO DEFINE WAKE LINES 


CARD 20 


(KWPATCH) (KWSIDE) (KWLINE) (KWPANl) (KWPAN2) (INPUT) (NODEWS) 


0 3 

Terminates root 
vortex 


CARD 19 


CARD 20 


(IDENTW 

1 


IFLEXW 

0 


IDEFW) 

0 


(WNAME) 

Connecting Sheet 


Note type-1 wake 
for connecting sheet 

(KWPATCH) (KWSIDE) (KWLINE) (KWPANl) (KWPAN2) (INPUT) (NODEWS) 
1 3 0 0 0 2 0 


CARD SET 21/22: TO DEFINE WAKE LINES 


CARD 20 


(KWPATCH) (KWSIDE) (KWLINE) (KWPANl) (KWPAN2) (INPUT) (NODEWS) 
0 0 0 0 0 0 5 

Last Wake 
Input 


1 



This example shows the simplest form of rotor wake input and 
provides an initially prescribed skewed cylinder form for the 
wake tube generated by copying the first filament input around 
the edges of the patch. This is adequate for cases where no body 
is in close proximity to the rotor. For cases where the fuselage 
presents substantial interference, a more involved wake 
specification is required. An example of this is provided in 
Figure 17 and below. The technique is identical to that 
illustrated above for type-1 wakes varying in a spanwise 
direction. For the example, the rotor disc has been divided into 
16 azimuthal sections (columns). Although the analysis in 
program VSAERO can cope with an initially prescribed wake 
filament passing through a body, it is better from the point of 
view of numerical stability if all filaments are prescribed so 
as to pass over the outside. In the situation pictured in 
Figure 17, the fuselage is narrow relative to the disc panelling 
and only the three filaments from the front of the disc, from 
columns 8, 9 and 10, need to be bent to pass around the body. In 
the illustration, filament 1 is copied around the edge to fila- 
ment 1 i Filaments 8, 9 and 10 are input individually. Filament 
11 is a return to the original trajectory and this is copied 
around until the wake closes at filament 17. The data set for 
this looks as follows. It should be noted that the filament is 
associated with the trailing corner of the panel in questions 
(corner 2 in the example shown). 



(KWPATCH) (KWSIDE) (KWLINE) (KWPAN1) (KWPAN2) (INPUT) (NODEWS) 

CARD SET 20 1 2 0 1 6 20 

Attachment to Input as 

columns 1-6 X and Z 

CARD SET 21/22: WAKE FILAMENT INPUT AS X AND Z RELATIVE TO SHED- 

DING POINT FOR FILAMENTS 1 TO 6 AND TO BE COPIED 
TO 7. 


CARD SET 20 1 2 0 7 7 00 

Last section of Copies 

"constant" wake previous 

section 


CARD SET 20 1 2 0 8 82, 3 or 40 

Filament 8 Input as 
required 

CARD SET 21/22: WAKE FILAMENT INPUT AS X,Z OR X,Y OR X,Z AS 

NEEDED FOR FILAMENT 8. 
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ORIGINAL PAGB 11 
OF. POOR QUALITY 


^^^LINE 1 DEFINED 



INPUT WAKE LINES 


COPIED WAKE LINES 


NOTE: FOLLOW SAME PROCEDURE 

FOR ROOT WAKE LINES. 


LINE 17 (CLOSE 16) 
COPIED 



| SIDE 2 

( x i SEPARATION SIDE 

U 2^JMKEJ_INE 

SIDE 1 


TYPICAL TIP-EDGE 
SHEDDING PANEL 


Figure 17. Rotor Wake Input Schematic with Fuselage Present 
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CARD SET 29 1 


2 


0 


9 


9 2,3 or 4 0 





I 


l 


CARD 

SET 

21/22: 

FOR FILAMENT 9 





CARD 

SET 

20 

1 

2 0 

10 

10 

2,3 or 4 

0 

CARD 

SET 

21/22: 

FOR FILAMENT 10. 





CARD 

SET 

20 

1 

2 0 

11 

16 

2 

0 

CARD 

SET 

21/22: 

FOR FILAMENTS 11 
MENTS 1 THROUGH 

THROUGH 

7. 

16. 

SHOULD MATCH 

FILA 

CARD 

SET 

20 

0 

0 0 

0 

0 

0 

0 


Copies 16 to 17 
and closes outside 
of wake 


P 


Wake definition continues with the input along sides 3 and 4 
of the patch displacing the wake filament as appropriate to pass 
around obstructions. 


Since rotor wakes are type-4 wakes, the initial wake 
strength must be defined with CARD SET 23, specifying V INNER and 
^OUTER* This is updated internally as the calculation proceeds. 
Unlike conventional type-4 wake situations, there is no need to 
use NORVEL (on CARD 8) for rotor patches since the disc panel 
boundary conditions are set internally using local loadings sup- 
plied by the rotor blade element calculation. 

In a conventional VSAERO data deck, the wake data completes 
the input string required for the aerodynamics calculation. 
However, if any type-4 rotor patches have been called, the 
program expects to read the input data set for the rotor blade 
elmement calculation. This is described in the next section. 
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5.0 ROTOR BLADE ELEMENT MODEL INPUT DESCRIPTION 

When the presence of a rotor is signalled by the insertion 
of a type-4 patch, the program requires that the data be loaded 
to construct the blade element model. This data set identifies 
the body patch involved; contains the controls which limit print 
volume and iteration cycles; provides the description of the 
blade geometry and mass properties including twist, planform and 
airfoil section; and defines the flight conditions. 


5.1 Rotor Patch Identifiers (Cards R1 and R2) 

i, 

Since the blade element model is built on the framework 
provided by the panel model, with disc panels and blade segments 
being matched, the first data items identify the rotor patch. 

CARD SET R1 provides the patch number and orients the disc for i 

the blade element velocity components. The parameter roll is the 
relative angle of rotor, zero degrees for a main rotor and ninety * 

degrees for a tail rotor. The roll angle is used simply to 
resolve the induced velocity components generated in the body ' 

axis system into the correct relation for the blade element 
analysis. Card R2 provides a general title for the rotor calcu- 
lation. 


5.2 Output Print Controls (Card R3) ; 

t 

t 

As an aid in rotor performance diagnosis, different levels j 

of printout are available. Ac the most detailed level, all of 
the blade section geometric onset flow and loading data are t 

available at each radial and azimuthal station for every step in 
the iteration cycle. This includes in the first group, avail- 
able at every radius and azimuth: J 

blade section radius, span, chord, geometric pitch, angle of | 

attack, yaw angle, inflow angle, Mach number, velocity . 

components in rotor control axis system, induced velocity j 

components in body axis system, section lift coefficient, 
and section drag coefficient. ! 

The second group includes the loading data available every radius j 

and azimuth. These are: » 

blade loading (lift/unit span), H-force loading, total 
torque loading, drag torque loading, lift torque loading, 
segment lift, segment H-force, segment total torque, lift 
torque, and drag torque. 

The third group includes integrated blade data available at 
each azimuth location. This is made up of principally blade 
lift, H-force and torques and includes the blade flapping 
parameters, the flapping angle and the first and second flapping 
derivatives. 
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Rotor totals are printed for each iteration and are not 
selectable. 

The output data selection is completed with an option to 
print the airfoil data as input and as interpolated at inter- 
mediate stations. It is recommended that restraint be exercised 
in switching on the print options since large volumes of output 
are generated. Irrespective of the print option chosen, the plot 
file contains all of the data for the final iteration cycle. The 
printout default values are all OFF and a particular group of 
output must be selected by inputting a value of one for the 
parameter. Card R3 is arranged in columns of 10 as outlined 
below. 


ISPNT1 
Card 3 0 or 1 


ISPNT2 
0 or 1 


IBPNTl IPT1 IPT2 

0 or 1 0 or 1 0 or 1 


Group 1 


Group 2 


Group 3 


Airfoil Data 


] 


z. 


u 

£ 


k 


Blade section data 
at every radius and 
azimuth 


Blade totals 
at each 
azimuth 


5.3 Iteration Controls . (Card R4) 


The rotor model used in the delivered verison of the program 
is for a fully articulated rotor and the performance calculation 
may proceed in either of two modes. In the first, direct mode, 
the blade control angles are preset and the calculation goes one 
cycle. In the second mode the rotor forces and moments are 
requested and the controls adjusted through three embedded itera- 
tion loops to produce the desired levels. The parameters entered 
on Card R4 control this process. These parameters, with MCOUNT 
controlling rotor moments, LCOUNT controlling rotor lift, and 
NFLAP controlling blade flapping, and hence control axis orienta- 
tion, may be set individually or together using default values. 


\ 5.3.1 Blade Flappiag 

The innermost iteration is the blade flapping cycle 
which steps the blade around the azimuth calculating the loads 
based on tie local conditions and the blade response to condi- 
tions at the previous step. The blade is assumed to be fully 
articulated. The default value of the controlling parameter is 4 
(four full azimuthal cycles allowed to stabilize flapping). If 
the blade flapping cycle is inhibited, by setting NFLAP-1 the 
blade flapping parameters, the control axis angle and the fore, 
aft and lateral flapping inputs must be included on Card Rll. 
Even if the default flapping cycle is selected, NFLAP-0 (set to 4 
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internally), it speeds convergence if values of control and 
flapping angles are set on Card Rll. 



uTi 


5.3.2 Rotor Lift 

The rotor lift is modulated using the collective pitch 
control. A total of six iterations are permitted if the default, 
LCOUNT»0 (set to 6 internally) are used. A starting value of 
collective pitch must be entered on CARd Rll. The search for a 
converged value of collective pitch uses a quadratic fit through 
the successively updated calculated thrust values, comparing at 
each step with the target value. Blade flapping equilibrium is 
re-established after each change in collective pitch. 


i 

'i 


* 


I 

) 
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5.3.3 Rotor Moments 

The rotor pitching and rolling moment loop entered on 
Card R4 produces the default, three iterations. If a search is 
ordered, target values of pitching and rolling moment must be 
loaded on Card R13 or zeroes will be assumed. If MCOUNT is set 
equal to 1, values of lateral and fore and aft cyclic pitch must 
be input on Card R12. 


5.4 The Blade Elemen t Model Geometry (CARD SETS R6, R7 , R8, 
R9, RIO) 

In the program the basis for the blade geometry and break- 
down into radial sections comes from the panel model. Having 
identified the appropriate patch, the panel corner points on the 
first column of the patch are studied (by the program). At the 
same time the number of rows (blade segments) is set, NR, and the 
number and size of the azimuthal increments determined, NC, from 
the number of columns on the patch. The panel geometry is 
available at this stage as coordinates in the global axis system. 
To convert these into blade radii requires the input of the rotor 
center of rotation and this is supplied on Card R5. The panel 
edges become the radial boundaries between the blade sections and 
the disc panelling and blade segmentation correspond. These 
radii form the basis for the rest of the geometry input data. 
Figure 18 illustrates how this procedure is carried out. The 
first radial station, defined by the innermost panel edge, is 
assumed to be coincident with the flapping hinge for the articu- 
lated rotors modelled here. The rest of the blade geometry 
information is input at the blade radial stations corresponding 
to the panel edges. 

Blade geometry is defined with three parameters entered on 
CARD SETS R6 , R7 and R8. These, respectively, are the chord, the 
twist and the leading-edge sweep. Each card set contains two 
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basic parts. They are: first, a card containing an indicator, 
tncpt which signals if a constant value is to be read (INSET 1) 
S? S if' values arl to be read at the NR+1 defining radial stations 
(INSET-0). If a constant twist is selected, the total twist 
should be loaded defined in the conventional sense relative to 
the 0.75 radius. 

The rotor description is completed by entering the blade 
mass distribution, required by the flapping calculation on Card 
R9 and the number of blades on RIO. 


L. (Cards Rll, R12, R13, R14) 


Cards Rll through R14 provide the input which set up and 
control the rotor performance. 


5.5.1 


(Card Rll) 


Card Rll provides the information which orients t 
control axis system in space, sets up the initial 
applies a starting value of collective pitch and sets up the 
correct onset flow and rotor rotational speed. 

If the default, zero, value for the control axis angle, 
ALPHAC, is used here, the program requires that system drag (the 
negative of the rotor propulsive force required) be loaded in 
Card R13. Since the program calculates the blade coning from the 
input mass properties, entry of the blade flapping values, A1 and 
Bl, completes the rotor orientation input. If the t? 

constrained, NFLAP-1 on Card R4, appropriate values of A1PHAC, A1 
and Bl, must be input. Otherwise, the calculation may be started 
assuming zero flapping, but supplying realistic starting values 
speeds the flapping convergence. Input of unrealistic values of 
ALPHAC, Al and Bl can lead to failure to close the flapping loop 
and a subsequent crash. 

Blade speeds are set with the input of the advance ratio, 
the rotor rotation rate, OMEGA, and the hover tip Mach number. 


5.5.2 


(Card R12) 


As was noted above in Section 5.3.3 in the discussion 
of Card R4, the user has the option, with the P* ra ® e ^ er ' 
of trimming the rotor to target moments (MCOUNT-O, default) ?[ 
setting cyclic controls and taking whatever moments result 
(MCOUNT-1) . If MCOUNT is set equal to 1, values of a and b 
must be entered on Card R12 or values of zero will be assumed. 
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5 . 5.3 


(Cards R13 and R14) 


Rotor lift and moment targets and the aircraft drag are 
entered on Card R13. If either LCOUNT or MCOUNT on Card R4 are 
left to the default value or if any value other than 1 is input, 
then target values of lift or moment must be loaded on Card R13. 

The value of drag entered here is the drag area formed by 
dividing the actual drag by the dynamic pressure and is used as a 
guide in setting up the control axis angle if this is not input 
above on Card Rll. 


Card R13 is completed with 
where it is derived. 


the entry of VcLIMB areas 


The air density used in the reduction of the loads and 
moments to coefficient form is entered on Card 14. If desired, a 
rotor tip loss factor may be used in the calculation. This is 
also entered on Card R14. When a tip loss is used, the lift 
outboard of the radius ratio, r/r TI p entered on Card R14 is 
varied linearly to zero at the tip. The calculated drag is not 
affected. 


5 *5.4 Airfoil Data Sets. (Card R15, CARD SETS R16) 

Airfoil data is used by the program in the conventional 
manner with table look-up and interpolation for C L an< j c D as a 
function of local aerodynamic angle of attack and Mach number. 
The airfoil sets are keyed to a particular radius, specified at 
input, and during execution the program interpolates between the 
data sets appropriate to the nearest radial stations on either 
side of the blade segment radius. 

At input the parameter, NDSEC, on Card R15 indicates how 
many data sets are to be loaded. Data sets may be loaded or 
copied from sets loaded earlier in the input string. This is 
indicated by a parameter on the first card of set R16. With the 
parameter, ICOPY-1, a data set is read. With ICOPY»0 the data 
set is copied form the previously read set. The radius (abso- 
lute) at which the data applies is also entered on Card R16.1. 
The data set follows using the "standard" C81 format. Ref. 2. In 
this format data is read as a function of blade section angle of 
attack for a range of Mach number. Each set of coefficients is 
entered separately. The program interpolates to the local coef- 
ficient value at the required value of Mach number and angle of 
attack and then between data sets, loaded as a function of span 
location, to the correct value. 






r 6.0 INPUT DATA DECK BLOCKING AND VARIABLE LIST 

; 

\ 

I t The VSAERO data deck assembly was described in great detail 

in the user's document and will not be repeated here since the 
only change in set-up from the operational point of view is the 
addition of the rotor data block. This enters the run stream 
after the wake information has been loaded. 


6.1 input Summary 

The input is divided into the following parts: 

(i) BASIC INPUT 

General information* operating mode* onset 
flow* reference conditions, special options 


(ii) PATCH GEOMETRY 

Description of configuration surface in com- 
|r ponents* patches* sections* basic points* 

etc.* for panel generation 


( iii) WAKE INPUT 

Wake-grid-planes* type of wake* wake separa- 
tion line* initial streamwise geometry 


(r) ROTOR INPUT 

Geometry* iteration and print controls* 
control settings, force and moment targets* 
and airfoil data 


(iv) SURFACE STREAMLINE INPUT 

Location of starting point for each surface 
streamline 

(v) BOUNDARY LAYER INPUT 

Reynold's number* etc. 

(iv) OFF-BODY STREAMLINE INPUT 

Location of starting point and required up- 
stream/ downstream distances for each off- 
body streamline 


In the following description* the input variables are first 
listed in 6.2 for each of the above parts. Then* 6.3 gives a 
detailed description of the function of each input variable. 
This is followed in 6.4 by an input flow chart to help with the 
assembly of the input data file. Section numbering has been left 
common with the original VSAERO document. In the detailed des- 
cription and flow chart sections, only the rotor input is des- 
cribed. The user is referred to the Program VSAERO User's Guide 
for the other sections. 
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6.2 input Variable List 




Basic Input Sununaiy 


Card No. 

Variables 

Format 

1 

Text 

20A4 

2 

IPRI, IPRLEV, IPRESS, MSTOP, MSTART , MODIFY 

615 

2A 

I PRGOM , IPRNAB , IPRWAK, IPRCPV, IPRPPI (only 
if IPRLEV-5 on CARD 2) 

515 

3 

MODE, NPNMAX , NRBMAX, ITGSMX, IMERGE, NSUB, 
NSPMAX, NPCMAX 

815 

3A 

NROWB(I), I-l, 1 NRBMAX ! (only if NREMAX< 0 on 
CARD 3) 

1615 

4(a) 

or 

NWIT, NVPI , IBLTYP (if MODE-1 on CARD 3) 

315 

4(b) 

NT, NHC (if MODE -2 on CARD 3) 

215 

4A 

(only if NVPI>0 and IBLTYP-0 on CARD 4(a) 


(i) 

NPSETS 

15 

(ii) 

NPCHBL, NBCOL, (KOL(I), 1-1, NBCOL) 

(Number of 4A(ii) cards - NPSETS) 

If MS TART >0 and MODIFY=0| this is the end of 
the basic data on a restart run. 

1615 

5 

RSYM, RGPR, RNF, RFF, RCORE, SOLRES, TOL 

7F10.0 

6 

ALDEG, YAWDEG, RMACH, VMOD, COMFAC 

5F10.0 

6A 

ALBAR, RFRECU, HX, HY, HZ (only if MODE-2 on 
CARD 3) 

5F10.0 

7 

CBAR, SREF, SSPAN, RMPX, RMPY, RMPZ 

6F10.0 

8 

NORSET, NVORT, NPASUM, JETPAN, NBCHGE 

515 


3 




1 




Card No. 

Variables 

Format 

t 

8A 

(NORPCH(I), NORF(I), NORL(I), NOCF ( I ) / 
NOCL(I), VNORM(I), ADUB(I), 1-1, NORSET) 
(only if NORSETX) on CARD 8) 

515, 

2F10.0 

• 

8B ( i) VORT (only if NVORT>0 

on CARD 8) 

(ii) (RXV(I) , RYV (I)r RZV(I), 

I»l, NVORT+1) 

F10.0 

3F10.0 

1 

8C 

(NPSPCH ( I ) , NPSRF(I), NPSRL ( I ) , NPSCF(I), 
NPSCL(I), 1*1, NPASUM) (only if NPASUM>0 on 
CARD 8) 

515 


8D 

( JETPCH ( I ) , JETRF(I), JETRL(I), JETCF(I), 
JETCL(I), VIN(I), VOUT(I), 1*1, JETPAN) 
(only if JETPAN >0 ON CARD 8) 

515 

2F10.0 


8E 

(KPAN(I), KSIDE(I), NEWNAB(I), NEWSID(I), 1*1, 
NBCHGE) (only if NBCHGE>0 on CARD 8) 

415 


Card No. 

Variable 

Eormat 

9 

CTX, CTY , CTZ , SCAL, THET (component card) 

5F10.0 

9A 

CPX, CPY, CPZ, CHX, CHY, CHZ (only if SCAL<0 
on CARD 9) 

6F10.0 

10 

IDENT, MAKE, KOMP, KLASS, PNAME (patch card) 
(Note: If MAKE-0, go directly to CARD 16) 

415, 

6A4 

11 

STX, STY, STZ, SCALE, ALF, THETA, INMODE, NODES 
NPS, INTS (section card) 

f 


II 




Card No. 

Variable 

Format 

12(a) 

BY, 3Z , X (INMODE-1) 


(b) 

BX, BZ, Y (INMODE-2) 

3F10.0 

(c) 

BX, BY Z (INMODE-3) 


(d) 

BX, BY, BZ (INMODE-4) 


(e) 

TC, INPUT (INMODE- 5 or 7) 

F10.0, 

(f) 

H, INPUT (INMODE-6 or 8) 

15 

(g) 

BX, RAD, THET (INMODE-12) 

3F10.0 

13 

XRB, NINT (after options 12(e) and 12(f) 

F10.0, 

15 

14 

NODEC, NPC, INTC, MOVE (use with CARD 12 and 

30X , 


and 13) 

415 

14A 

NPCH, NSEC, IB, LB (if NODEC<0 on CARD 14) 

415 

14B 

XPIV, YPIV , ZPIV, HX, HY, HZ, ROT (if MCVE-1 
on CARD 14) 

7F10.0 

15 

THETA2 , THETA1 (only if NODES <0 on CARD 11) 

2F10.0 

16 

NPC, INTC, KURV, NPTIP, NODES, NPS, NTS 

35X, 


(special tip patch) (only if MAKE-0 

315, 


on CARD 10) 

10X, 

315 

16A 

(S(I), Y ( I ) , Z(I), 1-1, NPTIP (only if KURV>1 
on CARD 16) 

3F10.0 


Wake Input Summary 


Card No . Variables Format 


17 

X (wake grid plane stations) 

F10.0 

13 

NODE, NPC, INTC, MARK 


30X, 

415 

19 

IDENTW, IFLEXW, IDEFW, 

WNAME (wake card) 

315, 5X 
6A4 

20 

KWPACH, KWSIDE, KWLINE, KWPAN1, KWPAN2 , 
INPUT, NODEWS, IDWC, IFLXL, DTHET 

915, 

F10.0 
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Card No. 


Variables 


21(a) 

SWPY, 

SWPZ, 

X 

(if 

INPUT-1 

on 

CARD 

20) 

21(b) 

SWPX, 

SWPZ, 

Y 

(if 

INPUT-2 

ON 

CARD 

20) 

21(c) 

SWPX, 

SWPY , 

Z 

(if 

INPUT-3 

on 

CAP.D 

20) 

21(d) 

SWPX, 

SWPY, 

SWPZ (if INPUT 1 

*4 

on CARD 20) 


NODEWC, NPC, INTC 

VIN, VOUT ... (if IDENTW-4 on CARD 19) 


Format 

3F10.0 


30X, 

315 

8F10.0 


Card No. 


Format 


R1 

IPRCH , ROLL 

110, 

F10.0 

R2 

TITLE 

80A1 

R3 

ISPNT1, ISPNT2 , IBPNT1 , IPT1 IPT2 

8110 

R4 

MCOUNT, LCOUNT, NFLAP 

3110 

R5 

XO, YO, ZO 

3F10.0 

R6.1 

INSET 

110 

R6.2 

CHORD (I), 1-1 for INSET— 1 
I-NR+1 for INSET-0 

8F10.0 

R7.1 

INSET 

no 

R7.2 

TWRATE , for INSET— 1 

TWIST (I) , I-NR+1 for INSET-0 

F10.0 

8F10.0 

R8.1 

INSET 

no 

R8.2 

SWEEP ( I ) , 1-1 for INSET— 1 
I-NR+1 for INSET-0 

F10.0 

8F10.0 

R9 

BMASS 

F10.0 

RIO 

NB 

no 


‘*7 


. 


MU, ALPHAC, OMEGA, MTIP, COLL, A1 , B1 
ALPHAS (Not used), A1S, BIS 
BW, DRAG, PM, RM, VCLIMB 
RHO, TIPLOS 


Format 

8F1G.0 


8F10.0 


8F10.0 

8F10.0 


NDSEC 


R16.1 ICOPY, YR 


R16.2 TITLE, (NMACH(I), N ALPHA (I) , 1-1,2) 

R16.3 MACH ( J) , J*l, NMACH 

R16.4 ALPHAI, (COI(I,J), J=l,9) 

R16.5 (COI(J), J»10, NALPHA) 

Repeat R16.4 AND R16.5 NMACH times 

Repeat to R16.3 twice COI(I,J)=C Lf T=1 , 
Cd' I * 2 

Repeat to R16.1 NOSEC times 


110 , 

F10.0 


8F10.0 

8F10.0 


8F10.0 


Repeat Rotor Card Set R1 - R16 for each type-4 patch. 


Card No . 


Format 


F, KP, NS (compulsory input if IBLTYP-1) F10.4, 

(Place one card no. 24 for each streamline) 215 


F, KP, NS (end of surface streamline data) 2F10.4, 

215 
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FTTflWU- , 


RNB, TRIPUP, TRIPOP, XPRINT, XSKIP 

(CARD 26 only present if NVPI>0 on CARD 4(a)) 


££Ltmai 

5F10.0 


I 




Card No . 
26 


Card No . 
27 


Off-Bodv Velocity Scan Tnput Summary 
Variables 

MOLD, MEET, NEAR, INCPRI, INCPRJ, INCPRX 
(Start of each scan box. Finish the set with 
a blank card) 

XO, YO, ZO, NP (if MOLD-1 on CARD 27) 


28 

XO, 

YO, 

ZO, 

29 

XI, 

Yl, 

Zl, 

29A 

(ALI (I) 

, 1 = 

30 

X2, 

Y2, 

Z2 

30A 

(AL2 (I) 

, I s 

31 

X3, 

Y3 , 

Z3 

31A 

(AL3 (I) 

, I' 

32 

XI, 

Yl, 

Zl 


(if 

MOLD- 2 

33 

X2 , 

Y2 , 

Z2 


(AL3 (I) , 1-1, lNP3l)(only if NP3<0 on CARD 31) 


34 


34A 

34B 


NAL, NTHETA, NRAD (if MOLD=2 on CARD 27) 

(ALKI), 1 = 1, INAL I ) (if NAL<0 on CARD 34) 

(ALTHET (I) , 1=1, I NTHETA I J ( if NTH ETA <0 on 
CARD 34) 


34C (ALRAD(I), 1=1, I NRAD I ) ( if NRAD<0 on CARD 34) 
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Format 

615 

3F10.0, 

15 

8F10.0 

3F10.0, 

15 

8F10.0, 

3F10.0, 

15 

8F10.0 

7F10.0 

5F10.0 

315 

8F10.0 

8F10.0 

8F10.0 
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fi 

1 

:! 

■1 


i 


Rotor In put Data Deck Description 


Note : All integers are right justified. 


Card 

Columns 

Variable 

Format 

R1 

1-10 

IPATCH 

no 



Rotor patch number 



11-20 

ROLL 

Rotor roll attitude 
=0.0 for main rotor 
= 90.0 for tail rotor 

F10.0 

R2 

1-80 

TITLE 

80A1 


R3 1-10 ISPNTl 110 

*1, prints blade radial varia- 
tion of blade section data 

» 0 , suppresses print 

11-20 ISPNT2 110 

=1, prints blade radial varia- 
tion of blade loads , etc. 

*0, suppresses print 


21-30 IBPNT1 110 

■1, Prints rotor forces and 

moment summary at each azimuth 

-0, suppresses print 

31-40 TPT1 110 


*1, prints input airfoil data 
*0, suppresses print 
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i 

I 

J 




I 




rwi*- 




i 






Columns 

Variable 



41-50 IPT2 


no 


■I# prints interpolated airfoil 
at each blade section 

= 0 , suppresses print 


R4 1-10 MCOUNT 110 

Controls cyclic pitch/rotor 
moment loop 

=0, sets default *3 

*1, runs to set values of Als and 
Bis 


n-20 Lamm no 

Controls collective pitch/rotor 

lift loop j 

I 

=0/ sets default =6; recommend . > 

default or 2.3 

• i 

• ( 

21-30 NFLAP 110 \ 

Controls blade flapping/tip path 
plane loop 

i 

*0 , sets default *4; recommend *4 
for main rotors, *1 for tail 

rotors i 


R5 1-30 iO*. UL. ZQ 3F10.0 

Rotor center in body axis system 
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I 

B 






card columns 

R6 (SET) 

6.1 1-10 


6.2 


R7 (SET) 
R7.1 


R7.2 


R8 (SET) 
R8.1 


R8 ,2 


R9 


RIO 


1-80 


1-10 


1-80 


1-10 


1-80 


variable 


Format 


RT.&nE CHORD 


Inset — 1, constant chord; read 
only 1 value on 6.2 


110 


Inset -0, variable chord; read 
(NR+1) values on 6.2 


Blade chord values, dimensions in 
feet 


8F10.0 


RT.ADB TWIST 


Inset — 1, constant twist rate; read 
total twist from 7.2 


inset -0, variable twist; read (NR+1) 
values on 7.2 


Blade twist values 

Inset —1, total twist in Col. 1-10 


Inset-0, section twist relative to 
0.75 radius, in degrees 


rt.ade LEA nTNG-Fnr.F. SWEEP 


Inset —1, constant sweep; 8.2 
blank (sweep » 0.0) 


Inset -0, variable sweep; read 
(NR+1) values from 8.2 


Sv jep values 

Inset —1, blank card 


Inset -0, leading-edge sweep 
positive aft in degrees 


1-10 BLADE MASS 


Distribution assumed constant; 
dimensional in slugs/ft. 


1-10 NUMBER QL B L AD ES. 


110 


8F10.0 


110 




8F10.0 


F10.0 


110 
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— 7 nrr ~ 


Columns 


Format 


BQTOR PARAMETERS 

1-10 ADVANCE RATIO 

11-20 CONTROL AXIS ORIENTATION 

=0 if drag value on Card 14 used 

* value if tip path tilt entered 
(degrees positive aft) 

21-30 OMEGA 

rotor rotational speed (radians/ 
sec.) 

31-40 HOVER TIP MACH NUMBER 

41-50 INPUT COLLECTIVE pitch (degrees) 

51-60 A1 

Input fore and aft flapping 
(degrees) 

61-70 fil 

Input lateral flapping (Al and 
B1 can be input =0.0) 


F10.0 
F10.0 

F10.0 

F10.0 
F10.0 

F10.0 
F10.0 i 


CYCLIC CONTROLS 

1-10 SHAFT AXIS ANGLE (not used) 

11-20 Als 

Fore and aft cyclic input (deg.) F10.0 

21-30 filS 


Lateral cyclic input (deg.) 


F10.0 



« 


£aid 

Columns 

Variable 

ESLtmaL 

R13 


ROTOR LOADS AND MOMENT VALUES 



1-10 

GROSS WEIGHT (lbs.) 

F10.0 


11-20 

DRAG 




(D/q) square feet 

F10.0 


21-30 

ROTOR PITCHING MOMENT TARGET 




(set MCOUNT-3 on Card 3) 
(ft. -lbs.) 

F10.0 


31-40 

ROTOR ROLLING MOMENT TARGET 




(set MCOUNT*3 on Card 3) 
(ft. -lbs. ) 

F10.0 


41-50 

V CLIMB (ft. /sec.) 


R14 


MISCELLANEOUS 



1-10 

AIR DENSITY (slugs/Ft. 3) 

F10.0 


11-20 

TIP LOSS FACTOR 




Linear tip loss applied outside 
input value, R/RTip 

Default a 0 . 0 , no tip loss 

F10.0 

R15 

1-10 

NDSEC 




Number of airfoil data sets 
following 

110 

R16 


ATRFOTL DATA SETS 




Repeat ND sec. times 


Note : 

Repeat R1 

through R16 for each rotor, type-4, patch 

used. 
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* 

f 


r 


kF.RnnvNAMTr sect ton data . These data tables are 
input in the standard format currently used in 
the Rotorcraft Flight Simulation Program* C-81* 
Ref. 2. 

Repeat NDSEC times. 


Column Variable Value DfiSCCiPtiflJl 

1-5 ICOPY(IDSEC) 0 Complete aerodynamic tables are read 

in for this section 

ID Sectional aerodynamic data is copied 

over from previously defined section 
ID 

Note ; If ICOPY (IDSEC) >0, the rest of CARD SET R16 is omitted for 
this section. 


Columns 

Variable 

Value 

1-30 

TITLE 

ANY 

31-32 

NMACH(l) 

3-18 

33-34 

NALPHA (1) 

3-99 

35-36 

NMACH (2) 

3-18 

37-38 

NALPHA (2) 

3-99 

39-40 

NMACH (3) 

3-18 

41-42 

NALPHA (3) 

3-99 


ANY Alphanumerical title of sets of tables 

-18 Number of Mach number entries in C, 
table 


C^ table 

Number of Mach number entries in C d 
table 


C. table 
a 


NUIuDet OI nacn uuiuuct cin-ixco a.** 

table. Not used. Set-0 


C table. Not used. Set-0 
m 
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Note; CARD SETS R16C through R16E are repeated as a group three 
times for C C d and C m ; that is, in the following des- 
criptions 

K -1 C* 

K ' 2 C d 

K * 3 C Not used in rotor/body performance 

m calculation at this time 



8-14 

MACH ( 1 ) 

arbitrary 

Lowest Mach number 

15-21 

MACH (2) 

N 

Next 

highest 

Mach 

number 

22-28 

• 

MACHO) 

• 

n 

• 

Next 

highest 

Mach 

number 

• 

• 

64-70 

• 

• 

MACHO) 

• 

• 

n 

Next 

Highest 

Mach 

number 


Note ; Additional card may be required with same format to input 
NMACH(K) values of Mach numbers. 



Columns 

Variable 

Value 

Desc 

ZXlL 

L* 

>tion 

1-7 

ALPHA (K) 

arbitrary 

Angle of attack 

, degrees 

8-14 

COI (K , 1 ) 

N 

Coefficient 

at 

MACH (1) 

15-20 

COI (K , 2) 

N 

Coefficient 

at 

MACH (2) 

22-28 

• 

COI (K , 3) 
• 

n 

• 

Coefficient 

at 

MACH (3) 

• 

• 

64-70 

• 

• 

COI (K , 9) 

• 

• 

n 

Coefficient 

at 

MACH (9) 
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T 


Columns Variable Value 


8-14 

64-70 


COI (K r 10) arbitrary Coefficient at MACH(10) 


COI (K, 18) " Coefficient at MACH(18) 

1. CARD R16E included only if NMACH(K)>9. 

2. CARDS R16D and R16E repeated NALPHA(K) times for each 
angle of attack. 







A ’ 
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APPENDIX: SAMPLE CASE 


Contains: 1. Sample Input 


2 . Sample Output 


Note: The case has been chosen to be roughly representative of 

the H-34 rotor as tested on the Ames Rotor Test Module. 

The test module has been generated as a body of revolu- 
tion (see Section 3.1.2). The rotor has been run in a 
requested thrust , fixed cyclic mode with the control axis 
tilt input. This sets the tip path plane angle and de- 
termines the propulsive force. The output has been some- 
what truncated but a sample of each section has been re- 
tained. Only output specific to the rotor portions of the 
program is retained. For body output examples the user is 
referred to Reference 1. 
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